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ABSTRACT 



We analyze line-of-sight atomic hydrogen (H I) line profiles of 31 nearby, low- 
mass galaxies selected from the Very Large Array - ACS Nearby Galaxy Survey 
Treasury (VLA- ANGST) and The H I Nearby Galaxy Survey (THINGS) to trace 
regions containing cold (T < 1400 K) H I from observations with a uniform linear 
scale of 200 pc beam -1 . Our galaxy sample spans four orders of magnitude in 
total H I mass and nine magnitudes in Mb. We fit single and multiple compo- 
nent functions to each spectrum to isolate the cold, neutral medium given by a 
low dispersion (<6 km s -1 ) component of the spectrum. Most H I spectra are 
adequately fit by a single Gaussian with a dispersion of 8-12 km s -1 . Cold H I is 
found in 23 of 27 (~85%) galaxies after a reduction of the sample size due to 
quality control cuts. The cold H I contributes ~20% of the total line-of-sight 
flux when found with warm H I. Spectra best fit by a single Gaussian, but dom- 
inated by cold H I emission (i.e., have velocity dispersions <6 km s" 1 ) are found 
primarily beyond the optical radius of the host galaxy. The cold H I is typically 
found in localized regions and is generally not coincident with the very highest 
surface density peaks of the global H I distribution (which are usually areas of 
recent star formation). We find a lower limit for the mass fraction of cold-to-total 
H I gas of only a few percent in each galaxy. 



Introduction 



Dwarf irregular (dlrr) galaxies in the nearby universe are laboratories for studying the 
most fundamental properties of gas evolution and star formation. Large, mult i- wavelength 
studies of these systems are only just beginning. However, recent, high resolution sur- 
veys of large, spir al galaxies in the ultrav iolet (e.g., NGS - iGil de Paz et al.l 120071 ). op- 



tical (e.g., NFGS - Hansen] |2000|; ANGST - iDalcanton et al.l 120091) . infrared (e.g., SINGS - 



Kennicutt et~aDl2003h . and radio (e.g., THINGS - IWalter et al~1l2008l : HERACLES - iLerov et al 



20091 ) have given us insight into where and when stars form in a galaxy. They h ave also pro 



vided clues to the gas conditions surrounding sites of current star formation (ILerov et al. 



20081 ; iBigiel et al.l 120081). The conclusion is that stars form i n regions rich in cold, dense 



molecular material (jKennicuttl Il998t iKennicutt fc Evansll2012l ). These studies mainly focus 
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on high metallicity systems, however. How these results translate to low metallicity envi- 
ronments has not been fully explored, mainly due to the fact the most commonly observed 
tracer of molecular material, CO, is notoriously difficult to observe at low m etallicity (e.g., 
Tavlor et al1ll998l : barone et allboool ; iLerov et al.lbood : ISchruba et aDl2012h . 



Nearby dlrr galaxies offer the opportunity to study the star formation process in low 
metallicity systems at high spatial resolution. Studies of t he star formation rates in dlrr 



galaxies defined from UV, 24/im , and/or Ha emission (e.g., ILerov et al.ll2008t iBigiel et al. 



20081 ; iRoychowdhury et al.ll2009l ) have found a non-linear correlation with the local atomic 



hydrogen (H I) distributions. However, the dependence of star formation on m olecular gas 



Lerov et al. 


(2008 


) and 


Biriel et al. 



(120081 ). The star formation rate surface density is shown to correlate linearly with the 
molecular gas surface density in these galaxies. These results can be understood if one 



assumes that efficient star formation requires molecular hydrogen (e.g.. lKrumholz et al.ll2009 
and references therein). 

Based on theoretical studies, the very character of the i nterstellar medium (ISM) is likely 
to change at low metallicity (e.g.. ISpaans fc Normanlll997l ). I Glover fc Mac Low! (120111 ) mod- 
eled the so called X-factor which measures the relationship between amount of detected CO 
emission and the abundance of molecular hydrogen. They found that the a mount of CO drops 



subst antially at low metallicities, consistent with prior predictions (e.g., iMaloney fc Black 



19881 ). This deficit of CO at low metallicities is supported by observa tions showing rela- 



tively low or no detections of CO emission in nearby dlrr galaxies (e.g., iTaylor et al.lll998 
Barone et al.ll2000t ILerov et al.ll2005t ISchruba et al.ll2012l ). Consequently, this lack of met- 
als also results in lower amounts of polycycl ic aromatic hydrocarbon (PAH) emission (e.g., 
Engelbracht et al.ll2005l ; I Jackson et al.l 120061 ) . The PAHs are important to the formation of 
the cold ISM since they are critical in the shielding of UV and soft X-ray photons. 

The long standing problem of studying the molecular component in the low metallicity 
dwarfs continues. We expect dlrr galaxies to contain molecular hydrogen since they are 
forming stars. However, since direct observations of the molecular gas responsible for forming 
stars in large samples of dlrr galaxies via CO are currently not feasible, we must find a 
different tracer of star-forming gas. 

One intriguing idea is to find the immediate precursors of the molecular gas. The 
assembly of star forming molecular cloud s is generally believed to require cold H I (e.g., 
Wolfire et al.l 120031 ; iKrumholz et al.l 120091 ) . In our galaxy, cold H I clouds have been ob - 



served to surround and even intermix with molecular c 
through studies of H I narrow self absorption (HINSA 



ouds (e.g. 



Krco fe Goldsmith! l2010h 



Li fc Goldsmit 



ing technique to find cold H I gas was pioneered by Young fc L0I 



n 



1996 



■2003r) . One promis- 
19971 ) in a sample of 
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nearby, star forming dlrr galaxies. These authors used high angular and spectral resolution 
H I data of two nearby dlrr galaxies to decompose line-of-sight spectra into narrow and broad 
Gaussian components. The narrow-lined gas (a ~ 4.5 km s -1 ) was found only in specific 
regions within the H I disk while the broad-lined gas (a ~ 10 km s -1 ) was found along every 
line-of-sight observed. The narrow Gaussian component was attributed to cold (T<1000 
K) H I gas while the b road lined gas was assumed to be warm, neutral hydrogen (T>5000 
K). lYoung et al.l ( 120031 ) studied three more dlrr galaxies, finding that they, too, contained 
evidence for cold H I gas. Othe r authors have used this technique to discover cold H I in a 
small number of other galaxies. Ide Blok fc Walter! ( 120061 ) investigated t he H I distribu t ion in 
NGC 6822 and also found regions rich with cold H I gas. Recent work by lBegum et al.l (120061 ) 
found evidence o f cold H I gas in six dlrr galaxies from the Faint Irregular Galaxies GMRT 
Survey (FIGGS; iBegum et al.l 120081 ) . To date, cold H I in emission has been discovered in 
10 nearby dlrr galaxies. 

Within the limited data available, different measurements of the cold H I do appear to 



give c omparable results. Two galaxies, DDO 210 and GR8, from the sa mple of lYoung et al. 



( 120031 ) overlapped with that of IBegum et al.l ( 120061 ). lYoung et al.l (120031 ) used the Very Large 
Array (VLA) at a li near scale of ~200 pc beam -1 with a velocity resolution of 1.3 km s -1 
while IBegum et al.l ( 120061 ) used the Giant Metrewave Radio Telescope (GMRT) at a linear 
resolution of ~300 pc beam -1 with a velocity resolution of 1.65 km s -1 . Despite the use of 
different facilities and different spatial/spectral resolutions, both studies find similar results 
for both the spatial distributions and velocity dispersions of the narrow component. This 
agreement between the two studies suggests that the measurement of cold H I is generally 
robust, and not extremely sensitive to the exact observing parameters. 

The technique of ident ifying a cold neutral medium has also been used in nearby, high 
metallicity spiral galaxies. iBraunl (119971 ) isolated cold H I gas using relatively low spectral 
resolution (6 km s -1 ) imaging from the VLA for 11 of the closest spiral galaxies. He combined 
spectra from discrete radial bins and found that these combined spectra all had a narrow 
Gaussian core (FWHM < 6 km s -1 ) superposed onto broad (FWHM ~ 30 km s -1 ) Lorentzian 
wings. He also found that the cold H I gas is filamentary and is found in clumps, preferentially 
in the spiral arms. Since the majori ty of star formation in spiral galaxies occurs in molecular 



clouds within the spiral arms (e.g., iGordon et al.l 120041 ). it is not surprising to find the bulk 



of the cold H I associated with these features. 

T he purpose of our s t udy is to build upon the p revious results of 



19971 ). lYoung et al.l (120031 ). IBegum et al.l (120061 ). and Ide Blok & Walter! (120061 ) in order to 



Young fc Lol (11996 



provide limits to the locations and amount of cold H I gas in a large sample of 31 nearby, 
low-mass galaxies using a common spatial resolution. We describe our galaxy sample in £J2] 
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and our methods of signal extraction in §3j Our results are described in §1] and we end with 
a summary of our conclusions in §|5j 



2. Galaxy Sample, Observations, and Data Reduction 

Galaxies for this work were taken from t wo major H I survey s of nearby galaxies: The 



H I Nearby Galaxy Survey (THINGS; AW0605; IWalter et al.ll2008h and the Very Large Array 
- Advanced Camera for Surveys Nearby Galaxy Survey Treasury (VLA-ANGST; AO0215; 
Ott et al. in preparation). Our sample consists of 24 galaxies from VLA-ANGST and 7 
dwarf galaxies in the M81 Group from THINGS. The VLA-ANGST observations all have 
a high spectral resolution of 0.65 - 2.6 km s -1 . The 7 galaxies from the THINGS sample 
have velocity resolutions of 1.3 - 2.6 km s -1 . The high velocity resolution is crucial for 
detecting narrow velocity components. Other galaxies from the THINGS sample were not 
used because they had velocity resolutions of 5 km s" 1 , too coarse for this type of analysis. 
Our final sample contains 31 galaxies with distances ranging from 1.3 - 5.3 Mpc (average 2.9 
Mpc) and H I masses between 4xl0 5 - lxlO 9 M Q . 

All observations were made using the Very Large Arrajfj] (VLA) observatory. We briefly 
outline the basic reductio n procedure h ere, b ut note that both surveys followed similar 



recipes, described fully in IWalter et al.l ( 120081 ) and Ott et al. (in preparation). Standard 
AIPSo processing of spectral line data was followed. Phase, bandpass, and flux corrections 
were made using typical VLA calibrators. Two sets of data cubes were produced: a natural 
weighted cube with a typical beam size of ~10" and a robust weighted cube with a beam 
size of ~7". Both sets of data have rms noise values of ~1 mJy beam -1 in a single line- 
free channel. We use the robust weighted cubes in this work, which allows us to use the 
smallest possible uniform linear scale. Further, robust weighting produces minimal side-lobes 
allowing for a clean-beam which better approximates the dirty beam str ucture than t hat of 



the more complex structure seen in the dirty-beam of natural weighting (lBriggslll995f ). The 
better dirty-beam approximation of the robust weighting is ideal for the extended emission 
observed in H I data. 

The main difference between the data for the two surveys is due to the inclusion of 
updated Expanded VLA (EVLA) antennas into the array for the VLA-ANGST survey (the 
EVLA has since been renamed the "Karl G. Janksy Very Large Array"). Unfortunately, 



1 The VLA telescope of the National Radio Astronomy Observatory is operated by Associated Universities, 
Inc. under a cooperative agreement with the National Science Foundation. 

2 The Astronomical Image Processing System (AIPS) has been developed by the NRAO. 
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the conversion of the digital EVLA signal to an analog signal (to be compatible with the 
VLA signal) aliased power into the first 0.5 MHz of the baseband. This only affected the 
EVLA- EVLA baselines and they were subsequently removed from the data (see Ott et al. in 
preparation for full details). Aliasing did not correlate into the mixed VLA- EVLA baselines 
and the VLA- VLA baselines were unaffected. To compensate for the loss of the EVLA-EVLA 
baselines, additional observation time was spent on each source. Tabled] gives an overview of 
the major observational properties of our galaxy sample. The table colum ns are: 1) galaxy 



name , 2,3) RA and DEC, 4) distance in Mpc (from TRGB measurements in lDalcanton et al. 



20091 unless otherwise noted), 5) total H I mass, 6) absolute I?-band magnitude, 7) the 
major-axis length in arcminutes at the ~25 mag arcsec -2 leve l (except for KKH 98 and KK 



230, which use the Holmberg system of ~26.5 mag arcsec -2 ; iKarachentsev et al.l 120041 ). 8) 
inclination ( IKarachentsev et al.l 120041 ). 9) the 200 pc beam size (see §3.l|) . 10) the spectral 
resolution, and 11) the rms of the noise in the line-free channels of the data cube. While 
some metallicity estimates exist for the most ma ssive galax i es in our sample, most do not 
have measured abundances. Using the relatio n in lLee et al.l (120061 ) which relates Mb to the 
oxygen abundance (see also lBerg et al.l 120121 ) . we expect a range of 7.0<12+log(O/H)<8.2 
with a median value of 7.7 for the low mass galaxies in our sample. 



3. Identifying Narrow Spectral Components 
3.1. Methodology 



We ana lyze l ine-of-sight spectra following methods similar to those outlined in lYoung &: Lo 



(1996. 


1997 


) and 


Young et al. 


(2003) 



kernel to a common linear scale of 200 pc beam -1 . To do this we computed the angular size 
equivalent of 200 pc for each galaxy using the distances given in Table [U We then used the 
CONVL task in AIPS which incorporates the initial beam size, beam position angle, and 
requested output beam size to compute the required convolving kernel. We note that our 
circular beam sizes can sample slightly larger linear scales depending on galaxy inclination. 
Spectra were then extracted through every l'/5 pixel (1" pixels for NGC 247 and NGC 3109). 
Extracting spectra from every pixel over-samples the beam, but for our purposes, does not 
influence the results. Also, the beam sizes reported in Table [Hand the typical galaxy rotation 
speeds of 10-20 km s -1 are small enough that kinematic broadening is minimal, regardless of 
galaxy inclination. Note that spectral broadening would work to hide a narrow component, 
not create one. 



Each extracted spectrum was successively fitted using our own routines with two dif- 
ferent functional forms: a single Gaussian and a double Gaussian. We also attempted to fit 
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the spectra with a fourth order Gauss-Hermite polynomial of the form 



[x 



ae 



h 



hi 



1 + -^=(2V2y 3 - 3V2y) + -£=(4y A - 12y 2 + 3) 



24 



where y = (x — b)/ (Jherm (Ivan der Marel fc Franxlll993l ). /13 and h± measure the amplitudes 
of an asymmetric and a symmetric deviation from an underlying Gaussian profile with am- 
plitude a, central velocity b, and standard deviation Oherm- G auss-Hermite polyn omials offer 
useful characterizations of non-Gaussian line profiles, but as [Young et al.l ( 120031 ) point out, 
the direct relationship between the Gauss-Hermite variables and the physical conditions of 
the gas is not obvious. In general, the results of the double Gaussian and Gauss-Hermite 
polynomial fits were in excellent agreement, although the detection efficiency of the Gauss- 
Hermite polynomials is much lower. However, we expect the obs erved spectra to have a 
Gaussian profile if the gas has an isothermal density distribution (jMerrifieldl Il993[ ) . Since 
the relationship between the Gauss-Hermite polynomial parameters and the physical nature 
of the gas is not obvious, we only report the results of the single and double Gaussian fits. 

A spectrum was identified as containing narrow-line H I only if multiple criteria were 
fulfilled. First, we only fit extract ed spectra wi t h a S /N greater than 10 (see §3.21) . We 
define S/N in the same manner as lYoung fc Lol (119961 ). that is, the peak in the spectrum 
divided by the rms noise in the line free channels of the data cube. We also required each 
individual component of the double Gaussi an fit to hav e a m in imum S/N o f 3.1. This 
minimum S/N requ irement was also used by lYoung & Lol (119961 ). lYoung &: Lol (Il997l ). and 
Young et al.l (120031 ) and is motivated by our desire to obtain significant detections. Also, 
each spectrum (single or multiple component) was required to have a velocity dispersion 
greater than the velocity resolution of the data (0.65, 1.3, or 2.6 km s -1 ) beyond the la 
errors in the fits. Requiring the velocity dispersion plus error to be greater than the velocity 
resolution ensures we are not fitting noise spikes in the data. Fortunately, none of the fitted 
spectra were affected by this criterion. Typical errors in the velocity dispersion of each 
single and double Gaussian component are roughly 0.7 km s -1 . We also do not enforce any 
restrictions on the central velocity offsets between individual components, leaving this as a 
free parameter. 

To quantify which function best fitted each spectrum, we first computed the variance 
of the residuals for each fit. The ratio of the variances between the single and double 
Gaussian fits were then compared statistically via a single-tailed F-test. A double Gaussian 
fit was determined to be statistically significant if the probability of improvement over a 
single Gaussian fit was 95% or greater. We were as conservative with our approach as 
possible in o rder to provi d e secure lower limits to the tota l amo unt of narrow-line H I in 
each galaxy. lYoung fc Lol (119961 . Il997l ) , lYoung et al.l (120031 ) , and iBegum et al.l (120061 ) were 
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less conservative in their approach, requiring only a 90% significance of improvement over a 
single Gaussian profile. Because of our relatively high S/N cut-off, the results of our fitting 
are relatively independent of the F-test criterion that we use, as long as it is above roughly 
70%. Since the number of detections of double Gaussian profiles increases very slowly as 
the cut-off is reduced from 95% down to 70%, we chose 95% confidence to establish a very 
firm lower limit on the total cold H I content. Establishing a realistic upper limit is more 
problematic and is highly dependent upon assumptions. 



We further adopt the criteria of Ide Blok fc Walter! (120061 ) that cold H I has a velocity 
dispersion of less than 6 km s _1 . A velocity dispersion of 6 km s _1 results in an upper limit 
to the gas temperature of ~1400 K (assuming (3/2)kT = (l/2)m p <7 2 ). Thus, we did not 
ascribe a second component as cold H I when the fits detected multiple components with 
similar velocity dispersions between 6-10 km s _1 or one component with a velocity dispersion 
of 6-10 km s _1 and another component with a much higher dispersion. The vast majority of 
spectra that were best fit by double Gaussian profiles with both velocity dispersions greater 
than 6 km s _1 were double peaked (in both low a nd high inclination gal axies) which is 
normally associated with expanding structures (e.g., [Brinks fc Bajajalll986l ). Furthermore, 



these higher velocity dispersions are not typical for the co 



with velocity disp e rsions of 2-6 km s 1 (Young &: Lo 



1996 



d H I described in the literature 



Braunl 119971 : lYoung fc Lolll997 



Young et al.ll2003t Ide Blok fc Walter! 120061; iBegum et al.l 120061); w e therefore rejected these 
locations as having narrow line emission. Ide Blok fc Walter! (120061 ) also noted that the bulk 
of the spectra within this category for NGC 6822 were double-peaked, leading them to make 
a cut at 6 km s^ 1 . Thus, having a cut at 6 km s -1 ensures that we are both reporting cold gas 
and eliminating a large fraction of the double peaked spectra. The double-peaked spectra 
eliminated by this criteria represent less than 1% of the total spectra and are not included 
in any of the following analyses. We also identify spectra that were sufficiently fit by a single 
Gaussian profile with a velocity dispersion of less than 6 km s _1 as cold H I. Limiting our 
cold H I identifications to gas with velocity dispersions less than 6 km s _1 ensures we are 
not misidentifying warm H I as cold H I. 

One possible scenario is that the broad component is the combination of multiple cold 
gas (narrow Gaussian) clouds. Given that the typical line-of-sight spectrum is best fit by a 
single Gaussian profile (see §1]), the multiple cold H I clouds would have to conspire in such 
a way to maintain this Gaussianity. While this scenario is physically possible, it requires a 
large number of cold HI clouds with distributions in velocities and amplitudes that would 
systematically combine to maintain a total Gaussian velocity dispersion equal to the observed 
velocity dispersion for the warm HI phase. High resolution absorption line studies would 
be required to distinguish whether the broad component is truly broad or a composite of 
many narrow components. However, if we assume that the galaxies can be approximated 
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by a thin disk we would expect that the line centroid of the cold H I to be at or very near 
the line centroid of the warm H I. Warm H I gas scale heights in dwarf galaxies tend to be 



a few hundred pc (e.g., IWarren et al.ll201ll ) compared to their few kpc diameters, so this is 
not a bad assumption. Offset velocity gas due to inclination effects would then contribute to 
the broadening of the warm (broad Gaussian) component and less so to the cold (narrow) 
component gas. Not surprisingly, the average central velocity offset between the components 
of the best fit double Gaussian profiles in our sample is 0.35 km s _1 , which is less than 
the velocity resolution of our data. Thus, the observations are most simply interpreted as 
a combination of cold gas near the central velocity (when observed) and nearly ubiquitous 
warm gas. We find it highly unlikely that the broad component is a combination of multiple 
narrow components. This leads us to the conclusion that the H I is dominated by the warm 
phase and that the cold phase represents a much smaller fraction of the ISM in our sample 
of dwarf galaxies. 

The top panel of Figure [1] shows an extracted spectrum (black) from the VLA- ANGST 
observations of Sextans A with single Gaussian profile (blue dashed), double Gaussian profile 
(solid red), and Gauss-Hermite polynomial fits (dotted magenta). The bottom panel shows 
the residual to the fit. This spectrum has a clear narrow peak that is not fit well by the single 
Gaussian profile. Both the double Gaussian profile and the Gauss-Hermite polynomial fit the 
spectrum statistically better than the single Gaussian profile at the 99.9% confidence level. 
Although the Gauss-Hermite polynomial is a better fit than the single Gaussian profile, the 
residuals clearly show that the double Gaussian profile is even better. 

In summary, we identify cold H I by the following process: 1) only spectra with a S/N 
greater than 10 were fit, 2) each spectrum was successively fitted with a single and double 
Gaussian function, 3) each individual component of the double Gaussian fit was required 
to have a minimum S/N of 3.1, 4) the velocity dispersion plus error had to be greater than 
the velocity resolution of the data, 5) double Gaussian fits were accepted if they improved 
the fit over a single Gaussian profile at the 95% or greater confidence level in a single tailed 
F-test, and 6) cold H I can be described by single or double Gaussian profiles with velocity 
dispersions less than 6 km s _1 . 



3.2. Simulations of Synthetic Spectra 

There are two factors which could lead to non-detections of existing cold H I. Cold 
H I could be missed simply because the column density is too low to reach our S/N require- 
ment, or it could be missed because a higher S/N is required to deblend the cold H I from the 
warm H I. As a test of our ability to extract multiple components from a given spectrum, we 
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produced three sets of 50,000 uniformly distributed synthetic spectra that spanned the S/N 
of our sample (10 < S/N < 60; one set each at 0.65, 1.3, and 2.6 km s _1 velocity resolution). 
Each spectrum contained two Gaussian profiles that are representative of the warm and cold 
H I gas; a "broad" Gaussian with abroad between 8 and 12 km s _1 and a "narrow" Gaussian 
with Onarrow between 3 and 6 km s -1 . We allowed for an offset of ±2.5 km s _1 between 
the line centroids. Velocity offsets larger than 2.5 km s _1 allow for much higher detection 
fractions since the spectra become more asymmetric. The observed cold H I line centroids 
are also very similar to the warm H I line centroids (as discussed in §3.ip . For each S/N, we 
randomized the amplitude of the cold H I Gaussian. We then chose the warm H I Gaussian 
such that the S/N was conserved. As for the analysis with our observations, each individual 
component was required to have a S/N > 3.1. A new set of Gaussians was produced if each 
component did not have a S/N > 3.1. We next added random noise to each spectrum and 
then passed each randomly generated spectrum through our fitting routine to gauge our 
recovery rate. 

The top row of Figure |2] shows histograms of the total number of synthetic spectra 
(solid black line) and the total number identified as requiring multiple components by a 
double Gaussian fit (solid red line) as a function of S/N. The red histograms also include 
the locations best fit by a single Gaussian with a velocity dispersion of less than 6 km s _1 . 
The bottom row shows the detection efficiency as a function of S /N for the double Gaussian 
fits (including the best fit single Gaussians with a velocity dispersion of less than 6 km s _1 ). 
Figure |2] shows that we never identify all double Gaussian profiles within the range of S/N 
values scrutinized and our recovery efficiency worsens as the velocity resolution decreases. 
The differences in recovery efficiency between the 0.65 and 1.3 km s _1 simulations are less 
pronounced than from 1.3 to 2.6 km s _1 . 

Table [2] shows the results of the simulations for each velocity resolution for the double 
Gaussian profiles that were identified. Column (1) is the velocity resolution, column (2) 
is the ratio between the input and extracted broad Gaussian amplitude (A 0}Sim / A^ extr ), 
columns (3) and (4) are the average differences between the simulated and extracted broad 
Gaussian central velocities (Av broad ) and velocity dispersions (Aa broad ), column is the ratio 
between the input and extracted narrow Gaussian amplitude (A n;Sim / A n)extr ), and columns 
(6) and (7) are the average differences between the simulated and extracted narrow Gaussian 
central velocities (Av nar ) and velocity dispersions (Aa nar ), respectively. These results show 
that when we do recover two Gaussian profiles, our routines accurately reproduce the input 
Gaussian parameters. 

In Figure [3] we plot the amplitude ratio of the broad and narrow Gaussian components 
(Abroad/ A narrow ) versus the ratio of the velocity dispersions ( abroad /oWrow) for the 1.3 km s 1 
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velocity resolution simulation. We show these ratios in four different S/N bins to understand 
where our fitting routines have trouble identifying the two Gaussian components. The grey 
points represent all of the simulated spectra, the red points are those spectra identified as 
containing two Gaussian components, and the blue points are the simulated spectra best fit 
by a single Gaussian with a velocity dispersion of less than 6 km s _1 (even though the input 
spectra contains both a broad and narrow Gaussian as described above). As is shown in Fig- 
ure [2j the recovery efficiency is dependent upon the S/N. The best fit single Gaussians are 
predominately in a regime where the narrow component dominates the spectra (amplitude 
ratios less than 1 and velocity dispersion ratios less than 2.5). We also have difficulty iden- 
tifying the narrow component when the broad component dominates the spectra in all S/N 
bins. Lastly, our routines do not pick up multiple components when the velocity dispersions 
are similar. These plots show some of the complexities in simulating and recovering multiple 
Gaussian components. When we do identify the spectra as containing multiple components, 
however, our routines accurately reproduce the input parameters. The simulations indicate 
that we can be confident of our detections in the observed data. However, the simulations 
also show that we can only reliably compute lower limits to the cold H I content. A complete 
census of the cold H I would require both higher S/N spectra and a detailed knowledge of the 
intrinsic distributions in Abroad /A narr0 w and Abroad /^narrow (in order to make incompleteness 
corrections). 

We also produced 50,000 synthetic spectra of single Gaussian profiles to quantify our 
false positives. These spectra were generated using a velocity dispersion between 8 and 12 
km s -1 with S/N values between 10 and 60. We used a velocity resolution of 2.6 km s^ 1 
since this gives the worst recovery results. With our acceptance criteria, we failed to find a 
single false identification of multiple components in our synthetic spectra. We do not start 
to see the potential for false detections until the confidence level is lowered to ~70%. These 
simulations give us confidence that our fitting results are robust to false detections. 



4. Results 

4.1. Cold H I Detections 

Each galaxy in our sample varies in terms of the total number of spectra and observed 
minimum, maximum, and average S/N (and column density). In total, we have analyzed 
roughly 4,100 independent lines-of-sight. Table |3] summarizes the observed spectral proper- 
ties of each galaxy. The columns indicate 1) galaxy name, 2) total number of independent 
lines-of-sight scrutinized (N t ), 3) average S/N of the extracted spectra (<S/N> t ), 4) the 
approximate column density at a S/N of 10 (Nm,min), 5) the peak column density (Nni^eofc), 
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and 6) the average column density (<Nhi>). All of the H I spectra in DDO 82, KDG 73, 
KK 230, and KKH 98 have S/N values below our threshold and as a result had no spectra 
analyzed, reducing our final sample to 27 galaxies. 

Table H] summarizes the results of our fitting routine. The columns are described as fol- 
lows: 1) galaxy name, 2) the average S/N of the spectra where cold H I is found (<S/N> C ), 
3) the areal filling factor of the cold H I defined as the ratio of the number of cold H I de- 
tections and the total number of scrutinized spectra {J-'jui), 4) the cold-to-total H I mass 
fraction defined as the ratio of the total summed column densities of the cold H I Gaussian 
profiles and the total summed column densities in the areas scrutinized (J^mass)^ 5) an es- 
timate of the upper limit to the cold-to-total H I mass fraction assuming each line-of-sight 
contains cold H I {J-^ ass ] see §4.51) . 6 & 7) the average velocity dispersion of the narrow and 
broad Gaussian components (<cr n > and <cr 6 >), and 7) the average velocity dispersion of 
the locations where a single Gaussian profile was sufficient to describe the spectrum (<er s >). 
Column 5 also includes all of the best fit single Gaussian profiles with a velocity dispersion 
of less than 6 km s -1 , while column 7 excludes them. 

We detected cold H I in 23 out of the 27 galaxies in our final sample. Figure H] shows 
the spatial distribution of the cold H I for each galaxy. The left panels show the total 
integrated H I intensity maps with contours of 10 20 and 10 21 cm -2 overlaid. The right panels 
have the H I surface density contours overlaid on the locations of the cold H I emission. 
The majority of the cold H I is clumped in localized regions. Typically, the cold H I is 
not spatially coincident with the very highest peaks in the total H I distribution in a given 
galaxy, although it is mainly concentrated in locations where the to tal H I column density 



exceeds the canonical threshold for star formation of ~10 21 cm -2 (jSkillman! 119871 ). Cold 
H I described in each of the previous studies also shows a preference for being located in 
regions of high column density, but not necessarily associated with the highest H I columns. 
The concentration of the cold H I in these regions is not unexpected given these location 
have the highest S/N values. The four galaxies in which we do not detect cold H I (DDO 99, 
MCG09-20-131, NGC 3741, and UGCA 292) show no distinguishing characteristics to give 
us insight as to why they are non-detections. Each of their total H I masses, M# values, 
distances, stellar disk sizes, and star formation rates (see §4.7p are similar to other galaxies 
in our sample that do show cold H I detections. 



4.2. Comparison with Previous Work 



Fortunately, our sample of 27 g alaxies had s o me o verlap with previous studies using 
similar methods to detect cold H I. lYoung et all j2003h detected cold H I in UGCA 292 
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and GR8 using the VLA while Begum et all (120061 ) detected it in DDO 53, MCG09-20-131, 
GR8 and M81 Dwarf A using the GMRT. These five galaxies are amongst the faintest and 
smallest galaxies in our sample and as a result, we detected cold H I in only 3 of these 5 
galaxies: DDO 53, GR8, and M81 Dwarf A. Our selection criteria failed to produce a cold 
H I signature in UGCA 292 or MCG09-20-131. 

The thr e e gala xies i n which we d i d dete ct cold H I are included in Figure HJ For GR8, 
Young et al.l (120031 ) and iBegum et al.l ( 120061 ) each found c old H I in both the northern and 
south-western portions of the galaxy while lYoung et al.l ( 120031 ) also found some evidence 
in the eastern portion of the galaxy. In contrast, we only found cold H I in the northern 
regio n of GR8. For M81 Dwarf A we detect cold H I only in the eastern side of the galaxy 
while IBegum et al.l ( 120061 ) find a hint of cold H I in the western side of the galaxy as well. 
For DDO 53 we found good overall agreement with the previous study, although our total 
detection area is smaller. 

The data quality between each survey is similar, yet they provide minor differences in 
their results. Our simulations suggest velocity resolution differences may play a role in the 
detection differences between the studies. The velocity resolutions in our sample for these 
galaxies are 0.65, 1.3, and 2.6 km s _1 for GR8, M81 Dwarf A, and DDO 53, respectively. The 
literature velocity resolutions are 1.3 and 1.65 km s _1 for GR8, 1.65 km s" 1 for M81 Dwarf A, 
and 1.65 km s _1 for DDO 53. The galaxy with the most similar cold H I distribution, DDO 53, 
has the worst velocity resolution in our sample, 0.85 km s -1 worse than the literature value. 
GR8 has been analyzed with three different velocity resolutions and each study finds slightly 
different results. We reduced our velocity resolution to 1.3 km s _1 for GR8 to compare to 
the other studies and obtained no difference in our results. M81 Dwarf A was analyzed with 
similar velocity resolutions, yet there still exists slight differences. Despite analyzing data 
from different observatories with different velocity resolutions, the general results are all 
similar. It seems likely that the differences in each gal axy arise from smal l diff erences in the 



select ion criteria and beam shape. For example, both lYoung et al.l (120031 ) and IBegum et al. 



(120061 ) used naturally weighted data cubes which, on average, produce higher S/N spectra. 
However, as noted in ^ the robust weighting for combined array data is ideal for studying 
extended emission due to the better clean beam approximation of the dirty beam in the 
CLEAN algorithm. The previous studies also used a cutoff of 90% in the F-test statistic and 
also allow the narrow Gaussian component to be larger than 6 km s -1 , finding values up to 
~8 km s _1 . 



-14- 



4.3. Are the Cold H I Non-Detections Significant? 

If each spectrum consisted of both a broad and narrow component and were not inher- 
ently a single Gaussian profile, we would expect to find far more locations with a narrow line 
signature than we actually do. The top panel of Figure [5] shows a histogram for every inde- 
pendent line-of-sight (black line), the expected number of identified narrow line detections 
(blue) (based upon our detection efficiency as calculated in §3.21 and the assumption that a 
narrow component exists at each location), and the actual number of narrow line detections 
(red) as a function of S/N. The bottom panel shows the detection fraction as a function 
of S/N. Figure demonstrates the large gap between the number of spectra for which we 
are sensitive to the presence of narrow H I and the number of spectra for which there are 
detections. Our expected recovery fraction is insensitive to small changes in the reasonable 
ranges used for the parameters. 

The lower S/N spectra have fewer fractional detections. As discussed in H'S.'2\ the lower 
fractional detections at lower S/N values are due, in part, to our recovery efficiency but may 
also be due to a mini mum total H I column density required for the appearance of cold H I. 



Kanekar et all (120111 ) recently observed H I in the Milky Way and determined a minimum 
column density threshold of 2xl0 20 cm -2 for the formation of a cold phase of H I, which 
is just at the column density limit of where our sample begins. The top row of Figure O 
shows histograms of all of the observed column densities in our galaxy sample (black), and 
the column densities where we detect cold H I (red). The bottom row shows our detection 
fraction function of column density. 

The small fraction of narrow line detections at high S/N compared to what would 
be expected clearly shows that the cold H I that can be identified with this technique is 
inconsistent with a ubiquitous distribution. It should also be noted that there could exist 
cold H I at every line-of-sight below a S/N of 3.1 that we would be insensitive to. However, 
it would be unphysical for the presence of gas to be related to our S/N criterion and not to 
the H I column density. Furthermor e, absorption line s t udies in the Milky w ay indicate that 



cold H I gas is not ubiquitous (e.g., iBegum et al.l 12010c iKanekar et al.ll201ll ). 



4.4. The Velocity Dispersions of Cold and Warm H I 

The typical velocity dispersion for the cold H I is ~4.5 km s _1 . The broad components 
and best-fit single Gaussian profiles vary by galaxy, but they generally have similar values. 
We expect the broad and best fit single components to be similar if they are both tracing 
the same warm H I. Figure [7] shows histograms of the narrow (red), broad (blue), and 
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single (black) velocity dispersions for each galaxy. The dotted vertical line denotes our cold 
H I cutoff limit of 6 km s _1 . The single component histograms have been scaled in order 
to show the narrow and broad component histograms more clearly. Generally, the velocity 
dispersions of the broad components overlap with those of the single components. However, 
the peak in the broad component is offset from the peak in the single component towards 
higher values. This behavior is also seen for the majority of the cases in the literature, and 
is probably due to the fact that it is easier to identify a narrow component when the broad 
component has a much larger value. 

In Figure we plot the average velocity dispersion, <a>, of the warm (black) and 
cold (red) H I gas as a function of radius for each of our different velocity resolutions. We 
have omitted the locations best fit by a single Gaussian with a velocity dispersion of less 
than 6 km s _1 . The average warm H I velocity dispersion declines with radius. Th i s tren d 



is similar to the total H I velocity dispersion versus radius seen by iTamburro et al.l ( 120091 ). 
The decrease in the velocity dispersion with increasing radius may indicate a decrease in 
turbulent energy supplied by the underlying stellar population as the radius increases away 
from the center of the stellar disk. 



4.5. The Areal and Mass Fractions of Cold F£ I 



Recent work on galaxy simulations of spiral disks have put limits on the predicted 



volum e filling factor of the different gas phases in the ISM (e.g., Ide Avillez fc Breitschwerdt 



20041 ). The cold (T < 1000 K) ISM occupies <20% of a galaxy's volume, depending upon 
the star formation (superno va) rate. The cold neutra l medium, however, may occupy only a 
few percent of the volume ( jMcKee fc Ostrikerlll977| ). The actual value will vary by galaxy 
and must also be sensitive to the available ISM coolants (e.g., C+, CO, dust, etc.) and the 
global gravitational potential. 

Without knowing the exact 3-dimensional structure of the galaxies in our sample, we 
cannot calculate the volume filling fractions. Instead, we calculate the areal filling fractions 
(Ffui), i.e., the ratio of the number of lines-of-sight with detected cold H I and the total 
number of lines-of-sight with S/N > 10 (see Table HJ column 3). We derive an average J r fui 
value of 9%. The true volume filling fractions based upon our detections are most likely 
even lower than these values since the scale height of the warm H I will be larger than that 
of the cold H I. The low filling fractions we compute are consistent with the ISM model of 
McKee fc Ostrikeri rtl977h . 



We also computed a lower limit to the total mass (flux) contribution of the cold H I for 
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eac h galaxy using the narrow co mponent of the Gaussian fits {Fmass'i see Table HI column 
4). lYoung fc Lol1l996l . Il997t ) and Ide Blok fc Walter! tood ) found that roughly 20% of the 
total H I mass is in the cold H I, which is significantly higher than the ~3% typical for 
our values of F^ss- ^ i s unclear as to the reason for this large discrepancy since the data 
are of similar quality. The biggest differences may lie in the higher S/N resulting from the 
natur al weighting these other studies use as compared to the robust weighting we employ. 
Also, Ide Blok fc Walter! use a different statistical test in defining the acceptance of cold 
H I detections. These authors accept a double Gaussian fit as statistically better than a 
single Gaussian fit if the \ 2 value is improved by 10%. Furthermore, our values are lower 
limits given our selection criteria and sensitivities. Figure [9]plots J-fui and J-^ ss as a function 
of Mb (left) and Mhi (right). No obvious trends arise between the areal filling fraction and 
mass fraction with absolute 5-band magnitude or total H I mass. 

If we assume (as noted in §4.3|) that each line of sight has cold H I just below our 
detection limits, we can estimate an upper limit to the amount of cold H I in each galaxy. 
To do this we have used a representative cold H I Gaussian of amplitude 3.1<r and velocity 
dispersion of 4.5 km s -1 at each location where we do not detect cold H I. We then add this 
"extra" mass to our detected cold H I mass to produce an upper limit to the cold H I mass 
fraction (J-'^ass'^ see column 5 in Table @J. These values of J-'mass can range up to 50%. 
However, since the cold H I is most frequently associated with higher values of Nhi (see 
§4.ip . it is unlikely that the true values of cold H I mass get as high as these upper limits. 

Another, perhaps more reasonable way to calculate an upper limit to the cold-to-total 
H I mass fraction is to find the typical fractional flux contributed to the cold component 
along each line-of-sight and assume this fraction holds throughout the sample. For locations 
where we find both warm and cold H I components, the cold H I typically is not the dominant 
phase in the ISM. The top row of Figure [TD] shows a normalized histogram of the flux ratio of 
cold and total H I defined as the ratio of the column density of the narrow Gaussian profile 
and the total column density along the line-of-sight. The cold H I constitutes ~20% of the 
total line flux for any given line-of-sight. The bottom panels of Figure [10] show this same 
ratio as a function of S/N for each velocity resolution. If we assume that every line-of-sight 
contains cold H I, then an upper li mit to the cold - to-to t al ma ss rati o for each galaxy is ^20% . 



This value is what is reported by [Young fc Lol (119961 . 119971 ) and Ide Blok &: Walter! (120061 ). 
However , this value is also l ikely too high since it requires the entire galaxy to contain cold 
H I gas. iDickey et al.l (120001 ) observed the Small Magellanic Cloud in 21 cm absorption and 
found that the mass fraction of cold H I to be less than 15% of the total H I content. It 
seems likely that the true value of the cold H I mass fractions of the low mass galaxies in 
our sample are between our lower limits and ~15%. 
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Begum et all (120101 ) observed 21-cm Milky Way absorption towards 12 background radio 



continuum sources. They were able to determine the flux contribution of the cold H I along 
many of these sight lines. They find that the cold H I contributes from a low of ~2% to as 
high as ~69%, of the total flux determined in the same manner, with a med ian of ^30% . 



While a direct comparison between the optically thick cold H I observed by iBegum et al. 



( 120101 ) and our optically thin emission is not obvious, we would not expect the fraction of 
cold H I in emission to be drastically different. Our result of ~20% is slightly lower than 
their median value, but does not appear to be at odds with their result. 



4.6. Cold H I Locations that Lack a Warm Component 

As discussed in §3.1[ some lines-of-sight were best fit by a single Gaussian profile with 
a velocity dispersion of less than 6 km s _1 . These best fit single Gaussian profiles were 
found in 6 of the 23 galaxies with cold H I detections. Most of these six galaxies have a 
large fraction of their H I at galactocentric radii beyond the ~25 mag arcsec -2 optical level, 
and the majority of these single Gaussian fits are also located at larger radii. Figure [11] 
shows the radial distribution of the cold H I locations that lack a warm component. The 
galaxies have been ordered from low (M81 Dwarf A) to high (NGC 247) absolute 5-band 
luminosity. Two of these ga laxies, M81 Dwarf A and Holmberg I, have large central holes in 



their H I distributions (see I Warren et all 1201 ll and references therein) where much of their 



gas is pushed beyond the optical radius. Only NGC 247 has a majority of these locations 
within the optical radius. 

Figure H2] shows an example of the difference between those areas that contain a warm 
component and those that lack one for NGC 4214. The left panel shows only those locations 
that contain both a warm and cold H I component while the right panel also includes 
the locations that lack a warm component. The contours correspond to the 10 20 and 10 21 
cm -2 total H I column densities and the blue circle approximates the 25 mag arcsec -2 
level (uncorrected for inclination and position angle). It is immediately apparent that the 
majority of the best fit single Gaussian p rofiles fall well beyond the bulk of the stellar 



population and far from any heat sources. iTamburro et all (120091 ) provides evidence that 



the warm H I velocity dispersions inside the optical radius of disk galaxies is driven by stellar 
processes, while areas beyond the optical radius may be driven by thermal broadening and/or 
magneto-rotational instabilities. These authors also point out that the velocity dispersions 
decline as a function of radius, falling below the typical 10 km s -1 beyond the optical radius. 
We also see this declining trend in the velocity dispersion with radius (see Figure [8]). These 
results are suggestive that cold H I can be the dominant constituent of the ISM only in the 
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lower UV radiation field beyond the ~25 mag arcsec 2 optical level. 



4.7. Comparing the Cold H I Gas Mass with Molecular Gas Masses and Star 

Formation Rates 



We next compare the cold H I gas masses to the molecular gas masses and star formation 
rates (SFRs) in each galaxy. To make these comparisons we compute the SFRs and estimate 
the amount of molecular hydrogen each galaxy contains. The relationship between the SFR 
and molecular hydrogen is discussed in detail in £JTJ We used the SFRs given by both FUV 
and Ho luminosities to estimate the amount of molecular hydrogen. 

Table [5] lists the SFRs derived from FUV and Ha luminosities taken from the liter- 
ature for each galaxy with a detected cold H I signature. The columns are: (1) galaxy 
name, (2) mass of the cold H I (M co idHi) computed using the lower limit mass fraction from 
Table H] multiplied by t he total H I mass given in Table (3) apparent FUV magnitude 
(m FUV ; iLee et al.ll201ll ), (4) the logarithm of the FU V luminosity (L eii v), (5) the FUV 



SFR (SFR^y) computed using the relation given in ISalim et all (120071 ). (6) the H 2 gas 
mass (M™ V ) computed from SFR^t/y assuming a 100 Myr timescale and a star formation 
efficiency of 10% {M^ v = SFR FUV x 100 Myr / 0.1), (7) the approximate fraction of the 
cold ISM attributed to the cold H I (M fuv) defined as M rn id Ni I (M rnldHJ + M^ v ), (8) 
the logarithm of the Ha l uminosity ( L w„\ iKennicutt et al.ll2008l ). (9) the SFR derived from 
Lho using the relation in iKennicuttl (119981 ). (10) the H 2 gas mass (M^i*) computed from 
SFR^a assuming a 5 Myr timescale and a star formation efficiency of 10% (M^ = SFRpuv 
x 5 Myr / 0.1), and (11) the approximate fraction of the cold ISM attributed to the cold 
H I {M.Ha) defined as M co i dHI /(M co i dHI + M^). Using our lower limit mass fraction from 
the spectra above S/N of 10 and extrapolating it over the entire galaxy provides a lower 
limit to the total cold H I mass. 

The expected relationship between the currently observed cold H I and the computed 
molecular gas masses is not immediately clear. The mass fractions of the cold ISM attributed 
to the cold H I is lower for H 2 masses computed from the longer timescale SFR^c/v (Table 
|5j column 7) than they are for those computed from SFR# Q (Table [5j column 11). The 
time scale to convert the cold H I into molecular hydrogen is likely dependent upon the 
metallicity of the environment (a longer time scale is expected for lower metallicities) and 
local UV radiation field (a shorter time scale is expected for lower UV radiation). Since dwarf 
galaxies have both low metallicities and different ranges of UV radiation strengths, it is not 
obvious which molecular gas mass estimate provides the best comparison. Unfortunately, 
our current data do not allow us to investigate these timescales further. 
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The top and bottom panels of Figure [13] show plots of the SFR^y (filled circles) and 
SFRhq (open squares) as a function of the cold H I mass and total H I mass, respectively. If 
the cold H I is related to the star formation in any direct way, we would expect to find a trend 
in this relation. Indeed, the higher the cold H I gas mass, the higher the SFR, although this 
same trend is seen with the total H I gas mass and may just be a consequence of the more 
H I gas available, the higher the SFR. The middle panel shows the SFR efficiency (SFE#j) 
defined as the SFR divided by the total H I gas mass as a function of the cold H I gas mass. 
The SFE/// is roughly constant over our range of cold H I masses. The difference s between 



the de rived FUV and Ha star formation rates has been previously reported in iLee et al. 



( 120091 ) and is due at least in part to the fact that the SFRs are averaged over different 
timescales. Environmental conditions may also influence SFR ff„ by leakage of the ionizing 
photons from the star forming regions (see lRelano et al.ll2012l for details). 



We can compare our molecular gas mass estimates de rived from the SFR s to a direct 
observational estimate of the molecular gas in NGC 4214 (IWalter et al.ll200ll ) to see if our 
values are reasonable. These authors observed CO emission in the inner ~l'6x 1'6 of the 
star forming disk of NGC 4214 (an area much smaller than the H I disk, see §4.61) and de- 
rived a molecular gas mass of 5.1 xlO 6 M & . This value is similar to the molecular gas mass 
estimate we get from SFR^q, of 6.1 xlO 6 M & . The molecular gas mass value computed by 
Walter et al.l was based on an assumed CO to H 2 conversion factor (the so called 'X-factor') 
derived from observations of Galactic environments (i.e., high metallicity). Recent work 
on th e X-factor has shown that the relation changes with decreasing mean visual extinction 



(e.g.. lLerov et al.ll2009l ; iGlover fc Mac Low!l201ll ; ISchruba et al.ll2012[ ). The lower metallicity 
environment of dlrr galaxies results in lower visual extinction values and thus a different X- 
factor value. Depending upon the true mean visual extinction value of NGC 4214, the actual 
molecular gas mass c ould be higher by a factor of a few to an order of magnitude or more. 
Schruba et al.l (120121 ) derive an X-factor yalue roughly 40 ti mes higher than the Galactic 
value for NGC 4214. Using this value, the IWalter et al.l (120011 ) molecular gas mass increases 
to ~2xl0 8 Mq, closer to our derived H 2 mass using SFRi?;yy of 1x10 s M . ISchruba et al. 



(120121 ) also derive a molecular gas mass of ~2xl0 8 M for NGC 4214. The factor of 2 dif- 
ference between the values most likely arise from our basic assumptions of a timescale of 100 
Myr for the FUV emission and/or the star formation efficiency of 10%. Unfortunately, the 
lack of CO detections for the few galaxies in our sample that have been obs erved for CO pro 



hibit us from repeating this exercise for other galaxies in our sample (e.g.. Ilsrael et al.lll995 



Verter fc Ho"dg3ll995l : lLerov et al.l bood : ISchruba et aDboid l Nevertheless, our estimates 
of molecular gas masses, although uncertain, are not unreasonable. Direct comparisons be- 
tween the CO detections in NGC 4214 and our cold H I detections is the focus of future 
work. 
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5. Conclusions 

We have observed line-of-sight H I emission spectra in 31 nearby, low metallicity dlrr 
galaxies in order to search for cold H I defined by a velocity dispersion of less than 6 km s" 1 . 
We have detected it in 23 of 27 galaxies after quality control cuts were applied. The cold 
H I may be the future sites of molecular cloud and star formation and to date, are the only 
way to potentially trace star forming gas at low metallicity. Based upon our observations, 
we find: 



The cold H I gas is found in localized regions and usually at tota l H I column de nsities 
above the canonical threshold of star formation of 10 21 cm -2 ( jSkillmanl 119871 ) . The 



cold H I in a given galaxy is also not typically associated with the very highest surface 
density gas (see Figures 0] & |6]) . 

The cold H I has a typical velocity dispersion of ~4.5 km s" 1 (T < 800 K). 

We derive an average value to the volume filling fraction of the cold H I of 9% (assuming 
our lower limit detections). 

We find lower limits to the cold H I ga s mass fractions of a few p ercent, consistent with 



some models of the multi-phase ISM ( iMcKee fc Ostrikerl 119771 ) . 



• The SFE#j is roughly constant as a function of the total cold H I mass over the 
observed gas mass ranges. 

• The cold H I contributes ~20% of the line-of-sight flux in locations where we detect 
both a cold and warm component. 

• Cold H I gas that lacks a warm component is typically found at radii larger than the 
25 mag arcsec -2 optical radius. 

Future work will investigate the relationship between the cold H I, local ISM conditions, 
and tracers of current star formation. 
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Table 1: Observed Galaxy Properties 
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° lDalcanton et all (120091 1: iKarachentsev et all (I2004f ); c MCG09-20-131 may have a greater 
distance than indicated due to ambiguities in measuring the TRGB. 
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Table 2: Simulated Spectra Results 
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Table 3: Individual Galaxy Properties 
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Table 4: Individual Galaxy Results 
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Table 5: Star Formation Rates and Molecular Gas Masses 
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0.91 


UGC 8833 


5.46 


16.64 


24.96 
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Fig. 1. — An example spectrum from Sextans A. The error on the points is shown in the 
upper right corner. This spectrum has a S/N of 37 and is clearly not fit well by a single 
Gaussian profile. The residuals to the fits are shown in the bottom panel. Both the double 
Gaussian profile and Gauss-Hermite polynomial statistically fit the observed spectrum better 
than the single Gaussian profile at the 99.9% confidence level in a single-tailed F-test. 
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Fig. 2. — Results from the synthetic spectra containing two Gaussian components randomly 
generated and processed with our fitting routine. Top Row: The total sample (black) along 
with the recovered profiles (red; includes best-fit single Gaussian profiles with a velocity 
dispersion less than 6 km s _1 for the three different velocity resolutions (0.65, 1.3, and 2.6 
km s _1 ). Bottom Row: Our recovery efficiency versus S/N for the double Gaussian profiles. 
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Fig. 3. — The ratio of the broad and narrow Gaussian amplitudes (Abroad/ A 

narrow J ^ 

function of the ratio of the broad and narrow velocity dispersions (obroad/cWrow) for the 1.3 
km s _1 velocity resolution simulation in four different S/N bins. The grey dots are all of the 
simulated spectra. The red dots are those spectra identified as containing two components 
by our fitting routines, and the blue dots are the simulated spectra that were identified 
only containing a single Gaussian component with a velocity dispersion of less than 6 km 
s _1 . The values in the parentheses are the recovery percentages for the best fit double (left) 
and single (right) Gaussians. We clearly identify more of the spectra at higher S/N values. 
Our routines have a harder time recovering the two Gaussian components when the velocity 
dispersions are similar. 
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Fig. 4. — Left: The total H I integrated intensity map. Right: The cold H I locations best 
described by a double Gaussian (red) and single Gaussian (blue). The contours for all panels 
represent the 10 20 and 10 21 cm -2 column density levels. The scale bar at the bottom right 
denotes a linear size of 1 kpc. The 200 pc beam is shown at the lower left. 



- 33 - 




-34- 








o 




(O'OOOSr) 330 



(o'ooozr) 33a 



(o'ooozr) 33a 







(o'ooosr) 33a 



(O'OOOJP) 330 



(o'ooozr) 33a 



Fig. 



IU — Continued. 
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Fig. 



IU — Continued. 
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Fig. 5. — Top: Histograms of the total observed sample for all galaxies (black), the expected 
number of cold H I detections given our detection efficiency and assuming every line-of-sight 
contained multiple components (blue), and the actual number of cold H I detections (red) 
as a function of S/N for each velocity resolution. The histogram bins were chosen such that 
each contains at least one independent line-of-sight. The large gap between the blue and 
red histograms demonstrates the cold H I identified with this technique is not ubiquitous. 
Bottom: The detection fraction as a function of S/N. 
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Fig. 6. — Top: Histograms of observed column densities of our entire sample (black) and 
where we find cold H I for each velocity resolution. Bottom: The detection fraction as a 
function of column density. Our detection fractions are relatively constant across all column 
densities observed. 
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Fig. 7. — Histograms of the velocity dispersions of the best fit single Gaussian profiles 
(black), narrow component (red), and broad component (blue). The single Gaussian profile 
distributions have been scaled down by the number in the parentheses in order to discern the 
other two components. The dashed vertical line denotes our narrow-line cutoff of 6 km s _1 . 
The velocity resolution (a) is listed below the galaxy names. Typically, the broad component 
velocity dispersions overlap with the single Gaussian profile velocity dispersions indicating 
they arise from the same gas phase. 
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Fig. 8. — The average velocity dispersion, <a>, as a function of radius for the warm (black) 
and cold (red) H I gas for each of our velocity resolutions. We have omitted the locations 
best fit by a single Gaussian with a velocity dispersion of less than 6 km s -1 . The averages 
are taken over 200 pc bins (the beam size). The error bars are the dispersions of the values 
in the bins. The velocity of the warm component decreases with radius indicating a possible 
decrease in turbulence as the radius increases from the main stellar body. Some of the 
errorbars at large radii for the 2.6 km s _1 velocity resolution overlap with values below our 
6 km s~ x cutoff. This is due to the dispersions being affected by a few larger values in these 
bins. 
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Fig. 9. — The areal fraction for the cold gas {J~fui) and the lower limit to the fraction of 
H I gas in the cold phase (J~!j£% ss ) as a function M# (left) and Mm (right). The areal filling 
fraction and mass fractions do not correlate with M# or Mm. 
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Fig. 10. — Top: Normalized histograms showing the contribution to the total line strength 
of the cold H I for each velocity resolution. Bottom: The average cold-to-total flux ratio as 
a function of S/N for each velocity resolution. The vertical error bars are the dispersions 
in the bins. The horizontal bar over each point shows the bin size from which the average 
was computed. For each panel we have omitted the locations where a single Gaussian profile 
with a velocity dispersion of less than 6 km s _1 best fit the data. The cold H I typically 
constitutes only 20% of the total line flux for locations which contain both cold and warm 
H I. 
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Fig. 11. — Radial distribution plots of the locations of the cold H I that lack a warm (<r > 6 
km s _1 ) component as a function of radius. The x-axis has been normalized by the 25 mag 
arcsec~ 2 radius (vertical dashed line). The galaxies have been ordered faint (M81 dwarf A) 
to bright (NGC 247) absolute 5-band magnitude. A significant fraction of the cold H I is 
outside of the optical radius of each galaxy, except for NGC 247. 
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Fig. 12. — Left: The cold H I distribution (red) of NGC 4214 for those locations best fit by 
double Gaussians profiles only. Right: The cold H I distribution (blue) for those locations 
best fit by single Gaussian profiles with velocity dispersions less than 6 km s _1 only. The 
contours represent the 10 20 and 10 21 cm -2 total H I column densities. The green circle 
approximates the 25 mag arcsec -2 optical level. The beam is shown at lower left. The 
majority of the cold H I detections described by a single Gaussian profile are outside of the 
bulk of the stellar distribution. 
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Fig. 13. — Top: SFR F( yy (filled circles) and SFR# a (open squares) as a function of the cold 
H I mass. Middle: The SFR efficiency (SFE HI ) defined as the SFR divided by the total 
H I gas mass as a function of the cold H I gas mass. Bottom: The SFR as a function of the 
total H I gas mass. The SFR shows a linear trend with both cold and total H I mass. 



